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Abstract 

Collisions of antiprotons with He-, Ne-, Ni-like, bare, and neutral uranium are studied theoreti- 
cally for scattering angles close to 180° and antiproton energies with the interval 100 eV - 10 keV. 
We investigate the Coulomb glory effect which is caused by a screening of the Coulomb potential of 
the nucleus and results in a prominent maximum of the differential cross section in the backward 
direction at some energies of the incident particle. We found that for larger numbers of electrons 
in the ion the effect becomes more pronounced and shifts to higher energies of the antiproton. On 
the other hand, a maximum of the differential cross section in the backward direction can also be 
found in the scattering of antiprotons on a bare uranium nucleus. The latter case can be regarded 
as a manifestation of the screening property of the vacuum-polarization potential in non-relativistic 
collisions of heavy particles. 

PACS numbers: 34.10. +x,34.90.+q,31.30.Jv,31.15.Ew 



I. INTRODUCTION 



The FAIR facility at GSI in Darmstadt will provide high-intensity antiproton beams at 
energies between 30 MeV and 300 keV at a magnetic storage ring and at energies between 
300 kev and 20 keV at an electrostatic storage ring. Further deceleration of antiprotons 
to ultra-low eV energies will be feasible via heavy ion trap facilities. This will enable a 
large variety of new experiments, including various atomic-collison type experiments with 
targets at storage rings. In particular, investigations of the antiproton scattering by a heavy 
ion at low energies are anticipated with the new GSI facilities. These investigations can 
give a unique possibility to observe an interesting phenomenon predicted in Refs. 
and termed the Coulomb glory. The phenomenon consists in a prominent maximum of the 
differential cross section (DCS) in the backward direction at a certain energy of the incident 
particle, provided the interaction with a target is represented by the Coulomb attraction 
of the nucleus (partly) screened by atomic electrons. Note, that the pure Rutherford cross 
section shows a smooth minimum at 180°. 

In the present paper we examine the Coulomb glory in collisions of antiprotons with He-, 
Ne-, Ni-like, bare, and neutral uranium (Z = 92). The ions are chosen to have fully occupied 
shells with n = 1 (U 90+ ), n = 2 (U 82+ ), and n = 3 (U 64+ ), respectively. The calculations 
have been performed using both semiclassical and quantum theory. Besides the screening 
potential due to atomic electrons, the vacuum polarization potential is taken into account. 
Atomic units {h = e = m e = 1) are used throughout the paper. 

II. THEORETICAL APPROACH 

We adopt the potential scattering model. This implies that neither excitation of atomic 
electrons nor ionization due to the collision is taken into account. The atomic electrons 
are regarded as a source of an electrostatic screening potential only. The use of the elastic 
scattering potential model can be justified for the problem under consideration since the 
energy of the incident antiproton (which corresponds the Coulomb glory effect) is quite low 
and not sufficient for the excitation of core electrons in the He-, Ne-, and Ni-like uranium 
ions. We also neglect the polarization of the atomic electrons by the antiproton. This 
polarization effect is small for highly charged uranium ions. For the neutral atom, both 
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polarization and inelastic processes can be significant. In the latter case, our treatment can 
be considered as a first approximation retaining the most important features of the Coulomb 
glory phenomenon. 

Two different parts are involved in our calculations. First, we need to construct the 
electrostatic potential of the target. This problem is solved with the aid of the density- 
functional theory (DFT) approach. For heavy atoms such as uranium, the relativistic version 
of DFT must be used. Second, we need to calculate the scattering phaseshifts from the target 
potential. This part is non-relativistic since the velocity of the antiproton is much smaller 
than the speed of light even in the vicinity of the nucleus. Moreover, the motion of the 
antiproton can be treated semiclassically, and the corresponding semiclassical methods can 
be used for the calculation of the phaseshifts. 



A. The scattering potential 

The effective potential V(r) experienced by the antiproton colliding with a heavy uranium 
ion can be represented as a sum of three local potentials: 

V(r) = V n (r) + V u (r) + V l j(r) J (1) 

where V n {r) is the potential of an extended nucleus, Vh(V) is the Hartree potential produced 
by atomic electrons, and Vy(r) is the Uehling potential taking into account the effect of 
vacuum polarization. The potential of an extended nucleus is given by 

Here p n is the nuclear charge density, normalized by the condition 

j d 3 rp n (r) = Z (3) 

to the nuclear charge number Z. We employ the Fermi model 

Pn{r) = — 77^ , (4) 

1 + exp[(r - r )/a\ 

where the parameter a is chosen to be 2.3fm/(41n3), while r and N are derived from the 
root-mean-square ouclea.cWe^sandaeno— on coition g ft 
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The interaction of the antiproton with the atomic electrons is described by the electro- 
static Hartree potential Vn(r) 



\r — r 
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(5) 

dr' r ri p(r') + 4n I dr'r'p(r'), 



r poo 



where p(r) is the total electron density, normalized via 

d 3 rp{r) = N (6) 



to the total number of the atomic electrons N. The electron density can be expressed in 
terms of the electron wave functions: 

b 

Here and fb are the upper and lower radial components of the relativistic one-electron wave 
functions in the shell b, and g& is the number of electrons in the shell b. The one-electron wave 
functions of heavy ions were obtained within the relativistic density functional theory using 
the local spin-density approximation with incorporation of the orbital-dependent Perdew- 
Zunger self-interaction correction The DFT calculations were performed employing the 
methods described in Refs. 0,0]. Finally, the Uehling potential is calculated according to 
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Vu(r) = -— I dr'r'p n {r') 
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x / dtii+i-^^n (8) 
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with a being the fine structure constant. 



B. Calculation of the phaseshifts and differential cross sections 

Since the expected kinetic energy of the antiproton is as low as a few hundreds of electron 
volts (the maximum velocity corresponding to the classical trajectory amounts to about 
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0.01c), the non-relativistic scattering theory can be applied. In the present paper, we make 
use of the partial wave expansion of the differential cross section 
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Here k is the momentum of the antiproton, v = —Z c mp/k is the Coulomb parameter (Z c is 
the charge of the core and in the case of He-like uranium Z c = 90), and Pi(cos6) are the 
Legendre polynomials. As one can see, the differential cross section (DCS) is a result of the 
interference between two contributions to the total scattering amplitude: the pure Coulomb 
(Rutherford) amplitude and the amplitude due to non-Coulomb (short-range) terms in the 
scattering potential. The phase shifts 6f are produced by the short-range part of the scat- 
tering potential. They can be expressed as a difference between the total phase shift Si 
corresponding to the angular momentum I and the Coulomb phase shift 6f 
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dr 



A way to calculate the phase shifts 8f without solving the radial Schrodinger equation is 

13|. Under this approach, the phase 



provided by the variable phase method [s[ [k], 11. \l2, 
shift is evaluated by solving a first-order differential equation. The method is robust and 
easy to implement. In our case, the differential equation to solve can be written as 

-2mpkv(r)r 2 

(12) 

x [cos Sf(k,r)F l (k,r) - sin Sf (k, r)G t (k, r)] 2 . 
Here 5f(k,r) is the variable phase which depends on r and v(r) is the short-range part of 
the scattering potential 

v (r) = V(r) + — . (13) 

Fi(k, r) and Gi(k, r) are the regular and irregular Coulomb wave functions, respectively fl^ |. 
Eq. (Tl2|) is solved with the initial condition 8f (k, 0) = 0. Then the long-distance limit of 
Sf(k,r) gives the value of the phase shift: 

Sf= hWf(fc,r). (14) 



For the energy of the antiproton as low as a few atomic units, the absolute value of the 
Coulomb parameter is large: \v\ ^> 1. That means, the motion of the antiproton can be 



treated in the framework of the quasiclassical approximation 151 ] . In this limit the Coulomb 



wave functions in Eq. (fT21) can be replaced by their asymptotic expressions through the 
Bessel functions [ijj]. An alternative way to calculate the quasiclassical phase shifts is to 
extract them from the quasiclassical radial wave functions. The phase shifts 8f are then 
expressed as a difference of the two integrals representing the phases of the quasiclassical 
wave functions in the total scattering potential and in the Coulomb potential —Z c /r: 



I dr\2m p (E-V(r))-( l + 1 / 2 ^ 




V \ I 2 



(15) 



where r and r c are the classical turning points for the two motions, respectively. Note that 
both integrals in Eq. (Tl5!) diverge as R — > oo but their difference does not, and the phase 
shift is defined correctly by the right-hand side of Eq. ffT5|) . In the case Z c = (neutral atom) 
6f are the total phase shifts corresponding to the angular momentum I, and Eq. (Tl5l) still can 
be used to determine these phaseshifts. To apply Eq. ( Tl2|) for the neutral atom, the Coulomb 
wave functions Fi(k,r) and Gi(k,r) have to be replaced by the spherical Bessel functions 
ji(kr) and ni(kr). In this case v{r) represents the total (short-range) potential. It is also 
worth to mention that the nuclear radius is much smaller than the distance corresponding 
to any of the classical turning points r or r c and, therefore, the annihilation probability 
during the collision is negligible. 



The quasiclassical approximation implies that the quantum number I is large 15]. Thus 
the phase shifts corresponding to small values of / may not be calculated accurately by Eq. 
(fT5]) (as opposed to Eq. ffl2l) ). However, this is not crucial since a large number of partial 
waves makes comparable contributions to the scattering amplitude (in our case - hundreds 
and even thousands), and the result does not depend significantly on a few partial waves 
with small I. We have performed calculations based on both Eqs. (Tl2l) and ( Tl5l) and found 
that the results agree with each other very well. 

The more partial waves in Eq. (J9j) contribute constructively in the backward direction, 
the larger is the maximum of DCS at 9 = 180°. On the other hand, the more "classical" 
is the motion of the incident particle, the more partial waves make important presence in 
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the wave function and the scattering amplitude. For the same energy, this is generally the 
case for heavier particles. Thus the experiments with heavy particles (antiprotons) could 
be more decisive in detecting the Coulomb glory than the experiments with light particles 
(electrons). 

III. NUMERICAL RESULTS AND DISCUSSION 

We have computed differential cross sections of antiproton-ion collisions for various en- 
ergies of the antiproton and several electronic configurations of the ion. We tested both 
Eq. (112]) and Eq. f|T5|) to calculate the phase shifts and found that they give very close 
results for the antiproton energies in the interval between a few tens electron volts and 10 
keV, where we can expect the Coulomb glory in backward scattering. In Fig. [1] we present 
the cross sections for He-like uranium. To facilitate a comparison between the results at 
different energies, the differential cross sections have been scaled according to 



The scaled Rutherford cross section becomes independent of energy and the ion charge: 



and equal to unity at 9 = 180°. Then the value of da/dQ at 9 = 180° represents the ratio of 
the ion DCS and the corresponding Rutherford DCS and can serve as a quantitative measure 
of the Coulomb glory effect. As one can see from Fig. [U for all the energies used in the 
calculations (5 eV, 20 eV, 100 eV, 500 eV, and 1 keV) da/dQ as a function of the scattering 
angle 9 has a maximum at 9 = 180° that indicates the presence of the Coulomb glory. 
However, the strongest effect is observed at the energy E = 100 eV when the scaled DCS 
reaches the value 4.1. The width of the maximum at this energy is about 2°. For lower and 
higher energies the effect becomes less pronounced. In the vicinity of the main maximum 
there are interference oscillations as described by Eq. (Q. With increasing antiproton energy, 
the frequency of the oscillations becomes smaller and the main maximum broader. 

The scaled DCS for Ne-like and Ni-like uranium are presented in Figs. [2] and El respec- 
tively. With increasing number of electrons in the ion, the range of the antiproton energies, 
where a prominent DCS peak exists at 9 = 180°, increases too. For Ne-like uranium, the 
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FIG. 1: (Color online) Scaled differential cross section da/dTL, denned by Eq. (I16p . on He-like 
uranium for different energies of the incident antiproton: (a) E = 100 eV, (b) E = 20 eV, (c) 
E = 500 eV, (d) E = 5 eV, and (e) E = 1 keV. The dashed line represents the scaled Rutherford 
cross section. 

Coulomb glory is best observed at the energy 300 eV with the scaled DCS equal to 90 at 
9 = 180°. For Ni-like uranium, the corresponding energy is 2 keV, and the scaled DCS 
reaches the value 609 at 9 = 180°. While the maximum at 9 = 180° becomes higher with 
increasing the number of electrons and the energy of the incident antiproton, its width does 
not change significantly and constitutes about 1°. 

For the uranium ions, we can compare the differential cross section with the corresponding 
Rutherford cross section to determine the presence of the Coulomb glory, and the quanti- 
tative expression is given by the scaled DCS (fTBj) . For the neutral uranium atom we need 
to find another criterion since no Rutherford cross section can be defined in this case. As a 



8 




90 



& 



174 175 176 177 178 179 180 

Angle (deg.) 



FIG. 2: (Color online) Scaled differential cross section dd/dTL, denned by Eq. (I16j) . on Ne-like 
uranium for different energies of the incident antiproton: (a) E = 300 eV, (b) E = 500 eV, (c) 
E = 200 eV, (d) E = 100 eV, and (e) E = 1 keV. The dashed line represents the scaled Rutherford 
cross section. 

reference DCS, we use the differential cross section itself, averaged over the interval of angles 
between 165° and 170°, where no prominent minima or maxima occur. Defined in this way, 
the average DCS (da/dfl) represents a characteristic value at a particular energy and angle 
range not affected by the Coulomb glory. Then the scaled DCS is defined as 



In Fig. H] the scaled DCS is depicted for the neutral uranium target and several energies 
of the incident antiproton. The results confirm the trend already observed for uranium 
ions: with increasing number of electrons the Coulomb glory effect is shifted towards higher 
energies. For the neutral uranium, the strongest maximum at 9 = 180° corresponds to the 
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FIG. 3: (Color online) Scaled differential cross section da/dQ, denned by Eq. (|16p . on Ni-like 
uranium for different energies of the incident antiproton: (a) E = 2 keV, (b) E = 1.5 keV, (c) 
E = 1 keV, (d) E = 2.5 keV, and (e) E = 3 keV. The dashed line represents the scaled Rutherford 
cross section. 

energy 7 keV with the scaled DCS equal to 320. 

Observation of the Coulomb glory effect in collisions of antiprotons with bare uranium 
nuclei can be of particular interest. Certainly, one cannot expect a large deviation from 
the background Rutherford cross section because of very short range of the non-Coulomb 
interactions due to finite nucleus size and polarization of vacuum; the smaller the radius 
of the potential, the less number of phase shifts 6f make noticeable contributions to the 
differential cross section (Q. However, if such a deviation is detected, it becomes a direct 
evidence for a screening property of the vacuum-polarization potential in non-relativistic 
collisions of heavy particles. Note that the finite nucleus potential is extremely short-range 
and influences very few first phase shifts Sf only. Thus any significant deviation from the 



10 



350 

S 300 
o 

i 250 
o 
o 



200 



§ 150 



5h 

^3 100 h 
-a 

1 50 
&0 







1 1 1 


1 1 1 1 1 


1 i 


1 


— 






/ V 


— 




/ c 


/ — 








/ e 



174 



175 



176 177 178 
Angle (deg.) 



179 



180 



FIG. 4: (Color online) Scaled differential cross section da/dQ, denned by Eq. (|18h . on neutral 
uranium for different energies of the incident antiproton: (a) E = 7 keV, (b) E = 5 keV, (c) -E7 = 8 
keV, (d) E = 2 keV, and (e) E = 9 keV. 

Rutherford cross section at 9 = 180°, which results from constructive interference of con- 
tributions with different /, should be mainly attributed to the polarization of vacuum. In 
Fig. we show scaled DCS as defined by Eq. ffTB]) for the energies of the antiproton 100, 
400, and 800 eV. The maximum of the scaled DCS at 9 = 180° exists for all three energies 
while the largest deviation from the Rutherford cross section corresponds to the energy 400 
eV and amounts to about 7%. This is, however, a rather large value, if we compare its 
magnitude with a typical QED contribution to dynamical processes with heavy ions. For 
instance, the QED correction to the DCS for the radiative recombination of an electron by 
a bare uranium nucleus does not exceed the level of about 2% {4, 16] . 
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FIG. 5: (Color online) Scaled differential cross section da/dCt, defined by Eq. (|16p . on bare 
uranium nucleus for different energies of the incident antiproton: (a) E = 400 eV, (b) E = 800 eV, 
(c) E = 100 eV. The dashed line represents the scaled Rutherford cross section. 

IV. CONCLUSION 



In this work we have studied the backward scattering of low-energy antiprotons by highly 
charged and neutral uranium. We found that a maximum in the differential cross section 
at the scattering angle 9 = 180° exists in a wide range of energies of the incident particle. 
However, at some energies the effect is enhanced. Classical mechanics describes this phe- 



nomenon as a combination o 
termed the Coulomb glory 
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ory and rainbow scattering at particular energies; it was 



|2j. Our quantum-mechanical calculations showed that the 
Coulomb glory can be observed for the energies of the antiproton within the range 100 eV - 
7 keV, depending on the electronic configuration of the ion. In general holds, the larger the 
number of electrons, the higher the energy where the effect has its strongest manifestation. 
In the case of Coulomb glory, the differential cross section at 9 = 180° can be much larger 



12 



than the corresponding background cross section. Actually, the ratio of these two quantities 
ranges from 4 for the He-like uranium to 609 for the Ni-like uranium. We also investigated 
possible manifestation of the effect in collisions of antiprotons with bare uranium nuclei. In 
this case the scattering potential differs from the Coulomb potential due to finite-nuclear 
size and vacuum-polarization effects. Both interactions are of very short-range that prevents 
large DCS values in the backward direction. However, some deviation from the Rutherford 
cross section does exist and can be increased by appropriate tuning of the antiproton energy. 
If experimentally detected, this effect can be regarded as an interesting manifestation of the 
screening character of the vacuum-polarization potential. 
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